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A general analysis of the statistical behavior of the envelope of a fading 
signal V(t)e'* U) is presented in this paper. The statistics include the 
probability P(V ^ L) that the amplitude V(t) will fade below a specified 
signal level L; the expected number N(L) of fades of V(t) below L per unit 
time; and the average duration t(L) of fades below L. The model for the 
fading signal is a constant vector plus a random interfering vector which 
represents the resultant of all the received extraneous signals and noise. 
The theoretical results agree with three empirically observed power relation- 
ships obtained in deep fades of nondiversity signals: P(V ^ L) <* L 2 , 
N(L) oc L and t(L) cc L. The theoretical results are applicable to a wide 
class of fading problems. The analysis includes the previous works of Rice, 
Nakagami, Norton, Vogler, Mansfield, and Short as special cases. 

I. INTRODUCTION 

A general analysis of the statistical behavior of the envelope V{t)e'^ t) 
of a fading signal is presented in this paper. Our principal interests are 
the probability, P(V ^ L), that the amplitude V(t) will fade below a 
specified signal level L* ; the expected number, N(L), of fades per unit 
time below the specified level L; and the average duration, t{L), of 
fades below L. These statistics are all functions of the signal level L. 

The theory presented herein has been developed to complement and 
extend the empirical results developed by my colleagues at Bell Tele- 
phone Laboratories ~ from their extensive experimental experience. 
Published data of other workers have also been considered. 

The previous theoretical works on the statistics of a fading signal 
often assume a complex Gaussian model for the fading signal. The 
theoretical support for this assumption is that, by the central limit 
theorem, the real and the imaginary parts of the sum of a large number 
of independent interfering signals will be approximately Gaussian. 

* More precisely speaking, in a long time period containing a large number of 
fades, the distribution P(V ^ L) represents the expected fraction of this long time 
period that the signal amplitude V will fade below L. 
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For tropospheric radio links, this model seems to be satisfactory. 
However, for line-of-sight radio links, the results of a short pulse 
experiment 7 and the angle-of-arrival measurements 8,9,1 indicate that 
the number of interfering signals is usually fairly small. Ray tracing 
theory also indicates that for typical line-of-sight radio links the number 
of paths contributing to multipath propagation is unlikely to be large. 
Furthermore, the theoretical results of the complex Gaussian model 
do not agree well with the experimental data on the statistics of fading 
signals of line-of-sight radio links, especially for certain overwater paths 
with severe fading. 

In this paper, we do not impose the restrictive assumptions of the 
complex Gaussian model. Rather, we simply model the fading signal 
Ve'* as a constant vector plus an interfering random vector; i.e., 

Ve'* = 1 + Re'* - 1+ « + &, 0) 

where R, 6, a, and are the amplitude, phase, real part, and imaginary 
part respectively of the interfering vector. The interfering vector is 
described by the joint probability density function /(a, 0) and represents 
the resultant of all the received extraneous signals, echoes, rays, and 
noise. The analysis applies for R and either dependent or independent; 
6 uniformly or nonuniformly distributed; a and /3 either Gaussian or 
nonGaussian. Thus, the results of this analysis may be applied to a 
wide class of fading problems. 

This paper treats the problem in three parts: The first is concerned 
with the amplitude distribution of V. The second considers the number 
of fades N(L) and the average fade duration t(L). The final section 
investigates several special topics including ra-distributions, chi- 
distributions, the Rayleigh distribution, Rice distribution, log-normal 
distribution and the sum of n unit random vectors. 

Appendix A is a list of symbols and their definitions. 

II. SUMMARY OF RESULTS 

(i) In spite of great variations in fading environment and test 
conditions, the experimental data 1 "*'"" 14 ' on P(V ^ L), 

N(L), and t{L) of most nondiversity* fading signals obey the 
following three prevailing power laws of deep fades: 



* The diversity signal is the output signal of a diversity combining system with 
two or more input signals. A "nondiversity fading signal" is a fading signal that is 
not a diversity signal. 
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(in 



P(V ^ L) cc L 2 

N(L) « L 

l(L) cc L 



for small L. 



(2) 
(3) 
(4) 



The theoretical analysis shows that if the probability density 
function f(a, /3) of the resultant interfering vector, Re' e = a + j/3, 
is a smooth function which is neither singular nor zero at the 
deep fade point (a = — 1, /3 = 0), then the statistics P(V ^ L), 
N(L), and i(L) of deep fades follow the three prevailing power 
laws (2), (3), and (4). The easily satisfied condition, a> > 
/(— 1, 0) > 0, is sufficient to obtain these functional relationships. 
(it) The set of power laws (2), (3), and (4) apply for R and 6 either 
independent or dependent, 8 either uniformly or nonuniformly 
distributed, (or, a and /3 either Gaussian or nonGaussian, either 
independent or dependent) as long as /(a, /3) is smooth. 
If /(a, /3) is singular at (a = — 1, = 0), then the theory predicts 
that for small L 



P(V g L) cc L 2M 

JV(/,) cc Z, 2 "- 1 

l(L) cc L 



1 > H > I 



(5) 

(0) 
(7) 



The exceptional behavior (5) consistent with y. = 1/2 has been 
observed experimentally on certain overwater radio links with 
severe fading. In this case, the resultant interfering vector 
contains the strong water-reflected ray as a dominant compo- 
nent. Therefore, the probability density function j(a, /3) has a 
singularity at the position of the dominant component vector.* 
(iv) If f(a, /3) has a zero at (a = — 1, /3 = 0) or is negligibly small 
at (a = — 1, /3 = 0), the theory predicts that for small L: 



P(V ^ L) * L 2 " 
N(L) cc L 2 "" 1 
l(L) cc L 



n > 1 



(8) 

(9) 
(10) 



The composite fading signals of diversity combining systems 
obey the set of power laws (8), (9), and (10). For overland radio 



* The complex Gaussian model, which assumes that /(a, 0) is a two-dimensional 
normal density function, is unable to explain the exceptional behavior (5) and (6). 
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links, the parameter n is equal to the order of diversity. In our 
experiments, the nondiversity fading signal of a relatively short* 
radio link with path length 15.87 miles has also shown the 
exceptional behavior described by (8), (9), and (10). 
(v) The theoretical results (4), (7), and (10) indicate that the power 
law, l(L) * L, for the average fade duration is more universal 
than those of P(V ^ L) and N{L). This prediction agrees with 
available experimental data. 
(vi) In general, the relationship between f(a, /3) and the amplitude 
distribution P(V ^ L) is not unique. As an example, this non- 
uniqueness shows that specifying the Rayleigh distribution for 
the amplitude of a fading signal does not necessarily imply that 
a and /3 are Gaussian, nor does it necessarily imply a large 
number of interfering signals. 



* At 4-GHz operating frequency, the average path length of line-of-sight radio 
links is about 27 miles. 
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I. INTRODUCTION AND SUMMARY 

In the study of fading signals due to multipath interference, the 
experimental data on the cumulative amplitude distribution, P(V ^ L), 
of deep fades are often plotted on a graph paper where the fade depth 
is expressed in dB and the probability of fade is expressed on a log scale 
as shown in Fig. 1. The consensus based on large amounts of experi- 
mental data 1-0 '" -14,15-17 is that the cumulative amplitude distribution 
of most nondiversity fading signals in the deep-fade region can well be 
represented by a straight line with a prevailing inverse slope of 10 dB 
per decade of probability. 

The equation which describes this typical distribution on Fig. 1 is 

P(V ^ L) = e-L 2 , for L up ^ L ^ (11) 

where V is the envelope voltage of the random fading signal normalized 
to its nonfaded signal level, L is any specified signal level, e is a param- 
eter depending on fading environment, and L up is the upper bound of 
signal level below which the straight-line representation of P(V ^ L) on 
Fig. 1 is valid. 
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Fig. 1 — Cumulative amplitude distributions of fading signals. 



The empirical result (11) means that the amplitude distributions of 
most nondiversity fading signals obey the following square law of deep 
fades 

P(V ^ L) ex L 2 , L up ^ L ^ (12) 

in spite of the great variations of fading environment and test condition. 
However, there are some exceptional cases. The most profound excep- 
tion occurs on certain overwater radio links with severe fading. In 
these instances, the probability of fades, P(V g L), decreases very 
slowly as the signal level L decreases and is characterized by an inverse 
slope of 20 dB per decade of probability as shown in Fig. 1, implying 
a power law 

P(V ^ L) cc L (13) 

in the deep-fade region. 

Another kind of exception occurs on certain radio links with relatively 
little multipath fading. The probability of fades, P(V ^ L), decreases 
very rapidly with L and is characterized by an inverse slope of 5 dB 
per decade of probability as shown in Fig. 1. This kind of distribution 
follows the power law 

P(V ^ L) « V (14) 
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in the deep-fade region. For example, at 4-GHz operating frequency, 
the average path length of line-of-sight radio links is about 27 miles. 
The behavior (14) has been observed to occur on a relatively short path 
with path length of 15.87 miles. 

The theoretical amplitude distributions previously derived were based 
on a complex Gaussian model and predict a square-law dependence, 
P(V ^ L) cc L 2 , in the deep-fade region. For example, the Nakagami 
distributions, 18 which include Hoyt distribution, 1819 Rice distribution, 18 
and Rayleigh distribution as special cases, all are square law in the 
deep-fade region. The explicit expression of Nakagami distribution can 
be found in Equation (4.6-28) of Reference 18. Figure 1 also includes 
Rayleigh distribution (dashed line) for comparison with the experimental 
data. 

The small number of interfering signals in line-of-sight radio links 
suggests that the assumption of a complex Gaussian model may be 
unjustified. One of the main objectives of this paper is to determine the 
weakest set of assumptions under which the square law (12) is obtained, 
and the condition for which the exceptional case such as (13) or (14) 
will occur. 

In summary: 

(i) The theoretical model for the fading signal, Ve% is a constant 
unit vector plus a resultant interfering vector 

TV* = 1 + Re iB = 1 + a + jjS. 

The resultant interfering vector Re' e = a + j/3, with joint 
probability density function j(a, /3), represents the sum of all 
the received extraneous signals, echoes, rays, and noise. 
(ii) An infinite fade (i.e., V = 0) occurs whenever a = — 1 and = 0. 
At this point R = 1 and 6 = t. Therefore, the behavior of 
/(a, /3) near the infinite fade point (a = — 1, j8 = 0) is closely 
related to the power law of amplitude distribution P(V ^ L) 
in the deep-fade region. 
(Hi) For most radio links, the interfering signals and noise may be 
considered random, so that the joint probability density function 
/(a, /3) of the resultant interfering vector is a smooth function 
near the infinite fade point (a = — 1, /3 = 0). The analysis 
shows that if /(a, /3) is a smooth function which is neither 
singular nor zero at (a = — 1, /3 = 0), then P(V S L) <* L 2 for 
small L. The simple condition w > /(— 1, 0) > is easily 
satisfied by most radio links. The validity of this square law does 
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not require that /(a, /3) be a normal density function. Therefore, 
the number of interfering signals does not have to be large. 
(iv) The analysis shows that if f(a, 0) is not smooth, but is singular 
at (a = — 1, jS = 0), then for small L, 

P{V ^ L) « L 2 ", 1 > n ^ i 

A physical example for this case is the overwater radio link 
where the water-reflected ray is almost as stable as the direct 
ray. The resultant interfering vector in this case contains the 
water-reflected ray which is "not very random." In other words, 
the joint probability density function /(a, /3) has a high peak 
at the position of the dominant, stable component, and may be 
considered singular at that point. 
(v) If f(a, /3) is zero or is negligibly small at the infinite fade point 
(a = — 1,0 = 0), then the analysis shows that for small L 

P(V <: L) * L 2 ", It > 1. 

A physical example for this case is the short radio link where the 
phase differences among the multipath propagations are all 
small. Then the value of /(— 1, 0) is negligibly small because 
the phase, 6, of the resultant interfering vector is generally 
small. Another example for this case is the composite signal of 
the output of a diversity combining system where the artificial 
active combining device serves to create a zero of /(a, 0) at 
(a = -1,0 - 0). 

II. FADING SIGNAL MODEL 

The received fading signal is modeled as a constant vector plus an 
interfering random vector as shown in Fig. 2. The latter represents the 
resultant of all the received extraneous signals, echoes, rays, and noise. 
The received fading signal normalized to the magnitude of the constant 
vector can be written as 

V(t)e mt) = 1 + R(t)e mt) , (15) 

where R(t) and 0(0 are the normalized magnitude and the phase of the 
interfering random vector respectively; V(t) and <f>(t) are the normalized 
magnitude and the phase of the received fading signal respectively. 
Let x(t) and y(t) be the real part and the imaginary part of the com- 
plex fading signal Ve 1 *, i.e., 

7(fl«'* (1) = x(t) + jy(t); (16) 
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and let a(t) and 0(0 be the real part and the imaginary part of the 
complex interfering random vector Re' , i.e., 

R(t)e mt) = a(t) + m). (17) 

V, <f>, R, 9, x, y, a and /3 are all real random variables. The normalized 
output power is 

V 2 = 1 + 2R cos 9 + R 2 = (1 + a) 2 + /3 2 . (18) 

The relative phase, 9(t), between the interfering vector and the constant 
vector can be taken to have values from to 2r because 6 and (9 ± 2nir) 
for any integer n are indistinguishable to the received signal at an 
operating frequency. 

A geometrical interpretation of equations (15) to (18) shows that 
deep fades (i.e., small V) occur when R and 9 are near the infinite fade 
point, (1, tt), in the (R, 9) plane or equivalently when a and /3 are near 
the infinite fade point, (—1, 0), in the (a, 0) plane. 

For line-of-sight radio links, notice that as far as the received signal 
V(t)e'* (t) is concerned, scintillation, atmospheric divergence, and earth 
bulge effects may also be replaced by a mathematically equivalent 
interfering signal R(t)e iHn which when combined with the constant 
vector gives the received fluctuating signal. 

Therefore, the model described by equations (15) to (18) includes 
many possible fading mechanisms which may occur individually or 
simultaneously on a fading environment. 

Going one step further, this model also includes the situations where 
no physical constant vector* exists because the mathematical decom- 
position (15) is applicable to any arbitrary fading signal. For such 
situations, the constant vector may represent the average signal level. 
The fluctuation of the signal is considered to be caused by an equivalent 
resultant interfering vector Re' . 

For radio links subjected to multipath interference, the number of 
incoming component waves is usually more than two. Thus the resultant 
interfering vector Re' consists of more than one echo and should not be 
interpreted as a simple physical echo. The main distinction is that the 
magnitude and the time delay of a physical echo are not functions of 
operating frequency whereas the magnitude R and the equivalent time 
delay, f „ = 9/u, of a resultant interfering vector consisting of more than 
one echo, are functions of operating frequency (i.e., are dispersive). 



* For beyond-the-horizon radio links, there is no direct radio path between the 
transmitter and the receiver. 
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Fig. 2 — Fading signal model. 



III. GENERAL FORMULATION OF AMPLITUDE DISTRIBUTION 

Equation (18) shows that the probability that the random signal V 
be faded below a specified signal level L is equal to the probability that 
a and jS fall within the circular region 



(1 + «) 2 + (? ^ L 2 



(19) 



in the (a, 0) plane as shown in Fig. 3. Let /(a, /3) be the joint probability 
function of a and /3. Then P(V ^ L) is the integral of /(a, /3) over the 
circular region (19); i.e., 



P(V£L)= / f(cc, ^) 



da dp. 



(20) 



The statistical behavior of the interfering vector, Re' 9 = a + j&, is 
sometimes described by the joint probability density function q(R, 6) 
of the magnitude R and the phase 6 of the interfering vector. A similar 
derivation in terms of R and 6 yields 



P(V ZL)= / q(R, 6) 



dddR, 



where 



9 L = cos 



1+R 2 - £■ 
2R 



(21) 



(22) 



Most of the following analysis is in terms of a and /3. An equivalent 
result in terms of R and 9 is given in Appendix B. 



V 2 =0 + a) 2 +/3 




Fig. 3 — Fading signal model on (a, /3) plane. 
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IV. PROBABILITY DENSITY FUNCTION OF INTERFERING VECTOR 

Equation (18) indicates that the infinite fade (i.e., V = 0) occurs 
when a = — 1,0 = 0. For most overland radio links, the interfering 
signals and the noise vary continuously in a random manner, so that 
the joint probability density function j(a, 0) of the resultant vector is 
a smooth function. On the other hand, there are some paths for which 
/(a, 0) may not be smooth, but singular. For example, for overwater 
radio links, the water-reflected ray near the grazing angle is comparable 
in magnitude and stability to the direct ray. Therefore, the resultant 
interfering vector contains a dominant component and the probability 
density function /(a, 0) has a sharp peak (i.e., singular) at this point. 
If the heights of the antennas and the path length are such that the 
average phase of the water-reflected ray is equal to it, then the singu- 
larity of /(a, 0) occurs at the infinite fade point (a = —1,0 = 0). 

These discussions suggest that in the general analysis of deep fades, 
one should consider not only the case with smooth /(a, 0) but also the 
case where f(a, 0) is singular at (a = - 1, = 0). A general probability 
density function /(a, 0), which is useful in our study, is* f 

/(«, 0) = [(1 + «)' + PT^-Hia, 0) (23) 

= V 2( "- l) -H( a , 0). (24) 

where H(a, 0) is an arbitrary smooth function. In the range 1 > n ^ 1/2, 
the density function /(a, 0) has a singularity of order 2 | /* — 1 | at 

(a = -1, = 0). 

On the other hand, in the range <x> > y. > 1, the density function 
has a zero of order 2(/* - 1) at (a = — 1, = 0). For convenience, we 
shall call the parameter, n, the smoothness index of /(a, 0). 

Since the possible singularity or zero of f(a, 0) at (a — — 1, = 0) 
is taken care of by the factor [(1 + a) 2 + 2 ]*" 1 , we shall assume that 

oo > #(-1,0) > 0. (25) 

Thus, the density function /(a, 0) given by equation (23) is neither 
singular nor zero at (a = — 1, = 0) if, and only if, n = 1 because the 
condition n = 1 implies /(a, 0) = H(a, 0), and vice versa. 

In equations (23) and (24), the smoothness index n can be either an 



* The reason for the use of the factor (m - 1) instead of a simple power index in 
equations (23) and (24) is for the convenience of notation in Section IV of Part 3 
when we investigate the m-distributions. 

t The cases where the singularity of f(a, /3) occurs at positions other than the 
infinite fade point (a = —1, (3 = 0) will not be analyzed in this paper. A brief 
discussion is included in Section VII of this part. 
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integer or a noninteger. The only restriction on n is that 

M ^ i (26) 

The reason for this constraint on /* is that in Section V of Part 2 we 
find that if ju < 1/2, then the expected number of fades N(L) approaches 
infinity as the fade depth L approaches zero. This seems to be non- 
physical. Therefore, we require that n ^ 1/2.* 

Since H(a, fi) is an arbitrary smooth function and n can range from 
1/2 to oo, then the probability density function f(a, 0), as given by 
equation (23), includes a large variety of fading environments. 

V. POWER SERIES REPRESENTATION OF AMPLITUDE DISTRIBUTION 

We shall assume that H (a, fi) is sufficiently smooth so that the two- 
dimensional Taylor series 20 expansion of H(a, fi) is applicable in the 
neighborhood of a = —1, and fi = 0. Several situations, where the 
Taylor series expansion of H (a, fi) is not applicable, will be discussed 
in Section IX and Appendix C. 

The Taylor series 20 expansion of H (a, /3) gives 

H(a, 0) = E [A E c»ff B _,.,(-i, 0)(1 + ay- r fiT 1 

= H(-l, 0) + H ll0 (— 1, 0)(1 + a) + #o.i(-l, 0)0 (27) 
+ ± [ff 2l „(-l, 0)(1 + af + 2ff...(-l f 0)(1 + a)fi + H .i(-1, 0)fi 2 ] 

(28) 



+ ••• , 








where 








H tt . 


{ 1 t\\ TJl-u fl\ 






'•* l > 0) - ecr- r dfi rH(pt > p) 


a = -\ 






c n nl 






T r\ (n — r)V 



(29) 



(30) 

Substituting (27) and (23) into (20) for P(V ^ L) and carrying out the 
integration (Appendix D) yields 

P(V S L) = £ d 2S+2 L 2S+2 > (31) 



* Notice that the unity total probability requires that the singularity of /(a, 0) 
be integrable (i.e., p > 0). The constraint n ^ § does not violate this condition. 
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- d 2 L 2 " + £ d 2S+2 L 2S+2 >, (32) 

S = l 

where 

_ jjfltiia (33) 

M 



^2a-2y.2>(~"l> Q) 



(34) 



a2S+2 ~ S\ (S + »)2 2S t=-o HOS-»)! ' 

5! = S(S - 1)0S - 2). . .3-2-1. 
The corresponding amplitude probability density function p(L) is 



S = oo 



p(L) = ^P(V £L) - E (25 + 2/0 4,+iL"^- 1 . (35) 

Oiv ,S = 

In the deep-fade region where L is small, the limiting forms of 
P(V ^ L) and p{L) are 

P(7 g L) > xg( ~ 1,0) L 2 " (36) 

£-•0 M 

and 

p(L) >2tH(-1,0)L 2 "- 1 , (37) 

£-•0 

where m ^ 1/2. The power law of deep fades for the three different 
cases (t) n - I, (ii) 1 > m ^ 1/2, and (tit) m > 1 with their physical 
fading environments will be discussed in the following sections (VI, VII, 
and VIII respectively). 

VI. PREVAILING SQUARE LAW OF DEEP FADES 

For the nondiversity fading signals of most radio links, j{a, 0) is 
neither singular nor zero at (a = — 1, /3 = 0). Then n = 1 and /(a, /3) = 
H(a, 0). Equations (31) to (37) under this situation become 



P(V g L) = Z d 2S + 2 L 2S+2 , (38) 
s=o 

= tt/(- 1, 0)L 2 + d,V + d 6 L« + • • • , (39) 

i(a, 0) = ff(«, 0), (40) 

d 2 = tt/(-1,0) = tH(-1,0), (41) 

, _ ^ V"» /2S-2».2>(~1) 0) /^O'S 

rf2S+2 - ( S +1)!2 2S & frOGS-,)! ' V ' 
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^-^(-1,0) = ^ f(a^) 



(43) 



P(L) = E (2-S + 2) d 2S+2 L 2S+1 , (44) 

P(V ^L) >tj(-1,0)L 2 , (45) 



L-*0 



and 



p(L) >2tt/(-1,0)L. (46) 



Equation (45) means that as long as /(a, /3) is neither singular nor zero 
at (a = — 1, |8 = 0), then the cumulative amplitude distribution in the 
deep-fade region always obeys the square law 

P(V <> L) « L 2 , L up ^ L^ 0. (12) 

Notice that this conclusion does not depend on any specific probability 
density function /(a, /3) for the interfering vector as long as f(a, /3) is 
smooth and oo > /( — 1, 0) > 0. The conclusion applies for a and /3 
either normal or not, either dependent* or independent, either with 
zero mean or with nonzero means. The magnitude R and the phase 6 of 
the interfering vector can be either dependent* or independent and 
can be either uniformly or nonuniformly distributed. Therefore, this 
conclusion covers a wide class of signal fading problems. 

Apparently the simple condition, oo > /(— 1, 0) > 0, is appropriate to 
the nondiversity fading signals of most radio links because the square law 
of deep fades is representative of the experimental data. 1_6 ' n_14,15_16t 

Notice that the first terms of equations (39) and (45), irL 2 , are the 
area of the two-dimensional region on (a, /3) plane bounded by the 
circle, L 2 = (1 + a) 2 + 2 , in which V ^ L, as shown in Fig. 3. 

The coefficient /(— 1, 0) in equation (45) has been observed to depend 
upon path length, operating frequency, path profile, and geographical 
factors. From the experimental data of a large number of radio links, 
it is possible to deduce an empirical formula of /(— 1, 0) as a function of 
these parameters. 22 ,23,24 ' 25 

In equation (12), the upper bound, L up , of signal level below which 



* This conclusion does not hold if the correlation coefficient between a and or 
between R and is unity because the joint probability density f(a, /3) becomes 
singular. 

t This theoretical result also explains an experimental fact that the observed 
amplitude distributions of atmospheric radio noise are also characterized by the 
square law (12) in the small amplitude region. 21 
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the square law applies also depends upon fading environment. Our 
experimental data show that L up of most line-of-sight microwave radio 
links is above 0.3 (i.e., above —10 dB). 

If /(— 1, 0) is negligibly small, so that the first term in the power 
series (39) can be neglected, then the second term, d A L*, dominates, 
with the result that the amplitude distribution follows the power law 
P(V ^ L) <* L*. We have observed this behavior on a short radio link 
in the signal range from -10 dB to -20 dB (i.e., 0.3 ^ L ^ 0.1). 
For fade depths deeper than -20 dB (i.e., 0.1 > L ^ 0) the quadratic 
term of (39) again dominates; and the transition region between 
P(V ^ L) « L 2 and P(V ^ L) <* L 4 occurs at about -20 dB for this 
short path. 

An obvious reason that /(— 1, 0) is small for short radio links is that 
the multipath length differences are mostly less than a half-wavelength. 

VII. A DOMINANT COMPONENT INTERFERING SIGNAL 

For this case, n is bounded by 1 > m = V 2 - % equation (36) 

P(V gl) > *H(-1,0) v „ i > M ^ i , (47) 

L-0 M 

and the corresponding power law of deep fades is 

P(V ^ L) cc L 2 ", 1 > M ^ h- (48) 

Since n is less than unity for this case, then as L — > 0, the probability of 
deep fades decreases more slowly than square law (12). Physically, this 
means that the deep-fade problem for these links is more severe. The 
experimental data of two oversea paths (shown as curve 2 in Fig. 16.5a 
and as curve 1 in Fig. 16.5b of Reference 15) follow this power law (48) 
of severe fading. 

It is true, however, that some overwater radio links still obey the 
square law (12) rather than (48). 15,17 The reason is that because of the 
geometry of the radio link, the singularity in the density function 
q(R, d) may occur, if it exists at all, at a position far away from the 
infinite fade point (R = 1, 6 = -k). Then in the neighborhood of the 
infinite fade point {R = 1, = x), the density function q(R, 9) (or 
equivalently /(a, 0)) may still be a smooth function. 

For overland paths, it is possible that an exceptionally calm and 
stratified atmosphere would also create a stable, dominant interfering 
signal over a sustained period. Then the joint probability density 
function f(a, 0) may also be singular at the position of this stable, 
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dominant interfering signal. Therefore, the results of this section on the 
power law of severe fading may also occur on an overland radio link. 

VIII. DIVERSITY SYSTEMS 

When n > 1, the power law of deep fades is 

P(V £L) > 7rg( ~ 1, ° ) L 2 ", h > 1, (49) 

P(V ^ L) « L 2 ", y. > 1. (50) 

As L decreases, the probability of deep fades decreases faster than those 
following the square law (12). Physically, this means the problem of 
fading for this case is less severe than those following the square law (12). 

The experimental data 1,2,5 of composite signals of the outputs of 
diversity combining systems show that the amplitude distributions 
of composite signals in the deep-fade region obey the power-law equa- 
tion (50) rather than the square-law equation (12). 

Since /i > 1 implies f(a, /3) has a zero at the deep-fade point (a = — 1, 
j3 = 0), these results show that the artificial active combining devices 
of diversity combining systems serve to create a zero at (a = ■— 1, /3 = 0) 
of order 2(n — 1) of the density function f(a, /3) of the equivalent 
interfering vector of the output composite signal. The value of n depends 
on the order of diversity. By comparing the power law (50) to the 
experimental data 1 ' 2,5 and the theoretical results on the diversity 
systems, we find that for most overland paths, the value of /i for the 
composite signal is equal to the order of diversity. 

IX. ONE-ECHO MODEL 

In the model described by equation (15) for the fading signal, if 
there is only one echo and if the magnitude of this echo is a constant, 
then R = A is a constant rather than a random variable. For this 
idealized case, the joint probability density function q(R, 6) of R and 6 
contains a delta function 

q(R, 6) = 8(R - A)W(d), (51) 

where W(6) is the probability density function of the random relative 
phase 8 between the echo Ae' B and the constant vector. 

For convenience, we shall call this specialized model the one-echo 
model. (In the literature, it is also known as the two-ray model.) 

In the analysis, we shall assume that A ^ 1. The case where A > 1 
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can be treated similarly simply by switching the roles of echo and the 
constant vector. 

Since the joint probability density function of the interfering vector 
of the one-echo model does not belong to the class of j(a, 0) discussed 
in Section IV of this part, the results of previous sections are not 
directly applicable. Nevertheless, substituting the density function (51) 
into the general formulation (21) for P(V ^ L) and integrating over R 
(see Fig. 4) yields 

P(V £L)= f * W(9)d6, (52) 

Jt-Bl 



where 



h = cos 



_, (I + A 2 - V 
2A 



(1 + A) ^ L ^ (1 - A). (53) 



Since the behavior of W(d) in the neighborhood of 6 = r is important 
for the analysis of deep fades, we shall assume that W{6) is smooth in 
this neighborhood so that the Taylor series expansion of W(6) is appli- 
cable. Then 



W(6) = W(t) + W 1 (ir)(6 - tt) + ^^ (6 - tt) 2 + 



(54) 



where 



wM = S w(e) 



n = 1, 2, 3, 



(55) 



Substituting equation (54) into equation (52) and carrying out the 
integration yields 



L*=1+2RC0S 0+R 2 . 




Fig. 4— The range of phase, (tt - B L ) ^ g (tt + L ), in which V ^ L. 
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(1 + 4)-£ L fe (1 - A). (56) 

If the magnitudes of the two vectors are equal, then A = 1, and equa- 
tion (56) becomes 



w^«-»Sow + 



2 ^ L ^ 0. 

(57) 



Since 



cos 



•■(>-£) —'[4 -a: 



cos 



1-f 



L 



(58) 



(59) 



Then the behavior of P(V ^ L) given by equation (57) in the deep-fade 
region is 

P(V ^ L) > 2W(t)L; (60) 



£->0 

.'. p(L) >2W(tt). 



(61) 



This result shows that as long as the probability density function W(0) 
of the random phase 6 is neither singular nor zero at = t, then the 
cumulative amplitude distribution of the one-echo model with R = 1 
always obeys the power law 

P(V £ L) <* L (62) 

in the deep-fade region no matter whether W(0) is uniform or not. 
Equation (57) shows that the nonuniform part of W(8) contributes 
only to the high-order terms of P(V ^ L) and does not affect the 
behavior of P(V ^ L) in the deep fade region. 

If the distribution of the random relative phase is uniform in (0, 2w), 
then 



W{6) = ± 

W„(t) = 0, n ^ 1. 
Equation (56) specialized to this case is 



P(V ^ L) = - cos" 

7T 



1 + A 2 - V 
2A 



)■ 



(63) 



(64) 
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If A = 1, then equation (64) becomes 

P(V = L) =^cos- 1 (l- y). (65) 

As far as the deep-fade region is concerned, equation (57) can also 
be written as 

P(V %L) = 2W(ir)L + 0(L 1+ "), v > 0; (66) 

.-. p(L) = 2W(t) + O(L'), (67) 

where O(L') is a symbol to denote the component which goes to zero 
at a rate equal to or faster than that of L n as L — > 0. 

Although the one-echo model of this section and the other two cases 
discussed in Appendix C do not exhaust all the situations where H(a, ft) 
is not analytic, the main objective is to show that the assumption of 
Taylor series expansion of H(a, /3) in Section V of this part is not 
strictly necessary for the derivation of the power law of deep fades. 

To unify the representations for all the cases considered in this paper, 
we shall rewrite equations (66) and (67) as 

P(V gL) = 2W(t)L 2 " + 0(L 2 " + <) (68) 

and 

p(L) = 2W(7r)L 2 "- 1 + 0(L 2 "- 1+ '), m = \, n>0. (69) 

Then the amplitude distributions of deep fades of all the cases discussed 
in Part 1 can be summarized as 

P(V ^ L) = d 2 L 2 " + 0(L 2 " + ") (70) 

and 

p(L) = 2n d 2 L 2 "- 1 + 0(L 2 "- 1+ "), /x ^ h v>0. (71) 

If H(a, @) is continuous at (a = — 1, = 0), then d 2 = vH( — l, 0)//i. 
If H (a, 0) is discontinuous at (a = —1,0 = 0), then d 2 = ttH( — 1, 0)/m- 
For a one-echo model, d 2 = 2W(ir) and n = 1/2. 
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Part 2. Expected Number of Fades and Average Fade Duration 
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I. INTRODUCTION AND SUMMARY 

In Part 1 we investigated the amplitude distribution of a fading 
signal. In a long time period, the cumulative amplitude distribution 
P(V ^ L) tells us the expected fraction of this time period that the 
signal will fade below any specified signal level L. However, P(V ^ L) 
does not tell us anything about the dynamic aspects of the fading 
signal. For example, a large number of short fades and a small number 
of long fades may have the same amplitude distribution. 

Some communication systems may tolerate the short fades but not 
the long fades. Furthermore, in the design of a diversity combining 
device, a distortion equalizer, or an automatic gain controlling device 
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to combat the fading problem, one needs information on the dynamic 
behavior of the fading signal. 

In Part 2 we present the results of our investigation on the expected 
number N(L) per unit time that the signal V(t) fades below a given 
signal level L; and the average duration t(L) of fades below L. 

The analysis is based on the general integral formulation of N(L) 
by Rice 26,27 and Vigants 4 and our results for P(V ^ L) in Part 1. Again, 
we do not impose the restrictive assumption of the complex Gaussian 
model so that the theoretical results may be applied to a wide class of 
fading problems. 

In the study of N(L) and t(L), the experimental data 1-6 for N(L) 
or t(L) are often plotted on a log scale as shown in Figs. 5 and 6. It is 
an experimental fact that the data for N(L) and t(L) can be well 
represented by straight lines on this kind of graph paper for fade depth 
deeper than —10 dB, as shown in Figs. 5 and 6. The slopes of these 
straight lines are directly related to the power laws of N(L) and t(L) 
in the deep-fade region. The experimental observations of N(L) and 
t(L) are summarized below: 

(i) The experimental data show that N(L) for most nondiversity 
fading signals obeys the power law, N(L) a L, in the deep-fade 
region. 
(it) For a short radio link from Villa Rica to Palmetto, Georgia, 
the N(L) of a nondiversity signal follows the cubic power law, 
N(L) ex L 3 , in the deep-fade region. 




-10 -20 -30 -40 

SIGNAL LEVEL £ = (20 LOG, L) IN dB 



Fig. 5 — Number of fades below signal level L. 
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SIGNAL LEVEL & = (20 LOG 10 L) IN dB 



-50 



Fig. 6 — Average duration of fades below signal level L. 



(Hi) The experimental data of N(L) of composite fading signals of 
most dual diversity systems also follows the cubic power law, 
N(L) « L 3 , in the deep-fade region. 

(iv) The available experimental data on average fade duration l(L) 
all obey the universal power law, t(L) °c L, in the deep-fade 
region. This includes the fading signals of nondiversity systems, 
diversity systems, long radio links, and short radio links. 

In summary: 

(i) Our theory indicates that if the joint probability density 
function f(a, 0) of the resultant interfering vector, Re' 9 = a + j@, 
is a smooth function which is neither singular nor zero at the 
infinite fade point (a = — 1, j8 = 0), then for small L 

N(L) ~ «yf(-l, 0)L « L, 

where a is a constant approximately equal to the average 
positive derivative of the amplitude of the fading signal in the 
deep-fade region. 
(ii) If /(a, 0) is singular or zero at (a = — 1, = 0) then 

N(L) ~ tH(-1, OKI, 2 " -1 ex L 2 "" 1 

for fi ^ 1/2 and small L. The cubic power law, N(L) <=c L 3 , of 
dual diversity systems and short radio links can be explained 
by this result when \i = 2, which means /(— 1, 0) is zero or is 
negligibly small. 
(Hi) The theory predicts that the average fade duration always obeys 
the power law t(L) <* L for small L no matter whether j(a, 0) is 
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smooth, singular or zero at (a = — 1, j8 = 0). This means the 
power law, t(L) <r L, is invariant with respect to variations of 
fading environment and diversity combinations of fading signals. 
This prediction agrees with the available experimental data. 
(iv) The theoretical results on N(L) and t(L) for the one-echo model 
are shown to be incompatible with the experimental data of 
most overland microwave radio links. Therefore, the one-echo 
model is not suitable for the study of the statistics of fading 
signals of these radio links. 
(v) For line-of-sight radio links at 4 GHz and 6 GHz, the average 
positive derivative of the amplitude of the fading signal is esti- 
mated to range from 2 X 10~ 3 to 4 X 10~ 3 times V tef per second, 
where V ie[ is the signal level when there is no interference. 

II. GENERAL FORMULATION FOR NUMBER OF FADES 

The general expression for the expected number of fades per unit time 
of a random signal V(t), below signal level L has been shown 26 ,27,4 to be 

N(L) = [ V " Vp(V, V) \ VmL dV, (72) 

J v-0 

where V = dV/dt, and p(V, V) is the joint probability density of 
V and V. For the sake of completeness, a brief derivation of (72) is 
included in Appendix E. 

The joint probability density function p(V, V) can be written in 
terms of conditional probability 28 as 

P(V, V) = Pl (V I V)p 2 (V), (73) 

where Pi(V \V) is the conditional probability density of V under the 
condition that the signal level is V; and p 2 (V) is the probability density 
of V. Substituting (73) into (72) yields 

N(L) = p 2 {L) [ V " V Pl (V | L) dV. (74) 

Let us define 

K(L) - 2 [. " V Pl (V | L) dV. (75) 

J v-0 

The physical meaning of the definition (75) is that 7+(L) is the condi- 
tional average positive derivative of V under the condition V = L. 
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The factor 2 in (75) is based upon the assumption that p^V \ L) is 
symmetric about V = 0. 

Since V(t) is a random fading signal, at a given signal level L, the 
value of its time derivative V is also random. In general, the conditional 
average positive derivative f+(L) is a function of signal level L. 

By using definition (75), equation (74) becomes 

N(L) = if+(L)p 2 (L). (76) 



Or, equivalently, 



N(L) = %V + (L)-^P(V ^L), (77) 



which indicates N(L) proportional to the conditional average positive 
derivative of the fading signal and to the probability density of fades 
at V - L. 

III. GENERAL FORMULATION FOR AVERAGE FADE DURATION 

In a long time interval, T, containing a large number of fades,* the 
expected total length of time that the random signal V(t) spends below 
a specified signal level L is 

t(L) = TP(V ^ L). (78) 

The expected number of fades below L in this interval T is TN(L). 
Therefore, the average duration of fades below L is 

t{L) TN(L) N(L) ™ 

Substituting (77) into (79) yields 

m = ^_ mv±D , (80) 

V+(L)j> p{v ^ L) 

Equation (80) shows that the average fade duration is inversely pro- 
portional to the conditional average positive derivative of the fading 
signal at V = L. 



* In our experiment on line-of-sight radio links, the typical time interval T is 
a whole summer of more than 100 days in which there are more than 500 fades 
below — 10 dB relative to the nonf aded signal level. 
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IV. ASSUMPTION ON CONDITIONAL AVERAGE POSITIVE DERIVATIVE 

Most existing theoretical work on N(L) assumes that V and V are 
independent so that 7+ (L) becomes a constant which is independent of 
signal level L. In this paper we include the situation for which V and V 
are dependent and assume that f + (L) can be expanded into a Taylor 
series in the deep-fade region (i.e., small L); then 

f + (L) = a + a x L + a 2 L 2 + a 3 L 3 + • ■ • , (81) 

where 

a = lim f + (L) = F~ + (0 + ), (82) 

£-0 + 

a, = lim ^y" + (L), (83) 

L-.0+ OLl 

1 3 2 — 

a 2 = —. lim -=s F+(L), etc. (84) 

The justification for this assumption is not trivial, and includes the 

following considerations: 

(i) The theoretical results based on this assumption agree with the 

available experimental data. 
(ii) For a complex Gaussian model, the conditional probability 
density function p t (V \ L) is known. Then with the help of the 
work of Rice 28,27 the integration indicated in (75) for t+(L) can 
be carried out in closed form. These explicit expressions are 
discussed in Section III of Part 3. The results of this model show 
that if the power spectrum of the Gaussian noise is symmetric 
with respect to the frequency of the sine wave (i.e., the constant 
unit vector), then V+(L) = a is a constant independent of L. 
On the other hand, if the power spectrum of the Gaussian jioise 
is asymmetric with respect to the signal frequency, then V+(L) 
is a function of L and the nonconstant terms in equation (81) 
cannot be omitted. 

The theoretical work of Clarke, 29 Ossanna, 30 and Gans 31 on 
mobile radio indicate that the power spectrum of the fading 
signal is generally asymmetric with respect to the received 
carrier frequency unless the straight line joining the base station 
and the mobile antenna is perpendicular to the velocity of the 
mobile and the antenna pattern is symmetric with respect to 
this line. Therefore, thework on asymmetric power spectrum, 
and hence nonconstant V + (L), is not purely academic. 



3236 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1971 



(Hi) It is known 18 that the correlation between any real random 
variable c(t) at instant t and its time derivative e(t + f) at 
instant (t + f) vanishes if f = 0. This is often used to support 
the assumption that V and V are independent and hence V + (L) 
is a constant. However, we know that the vanishing of correla- 
tion between V(t) and V(t ■+- f) at f = does not imply the 
independency of V and V unless V is normally distributed. 
Including high-order terms in equation (81) removes the assump- 
tion of independency of V and V and enlarges the applicable 
scope of this theory. 

(iv) In equations (82), (83), and (84) we define those coefficients of 
the Taylor series as the limits of V + (L) and its L-derivatives 
at L = + from the positive side. The reason is that V(t) is the 
absolute value of a fluctuating complex signal; i.e., 

V(t) = | x + jy | = | TV* | . 



When the complex fluctuating signal V(t)e i4,(n crosses zero, its 
absolute value V(t) may have a cusp at V = as shown in Fig. 7. 
Therefore, the derivatives of V(t) may not be well defined 
at V = 0. However, the limits of the derivatives at 7 = + 
from the positive side are well defined. 
(v) Fig. 7 also shows that although V = is a minimum of V(t), 




Fig. 7 — The amplitude 7(0 may have a cusp at V = even though the complex 
fading signal is a smooth time function. For convenience, the complex fading signal 
is plotted as a real function in this example. 
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V(t) = d/dt V(t) may not be zero at V = + . Therefore, in 
equation (81), the constant term a = l^ + (0 + ) does not vanish 
for most cases, and cannot be omitted from equation (81). 

V. POWER SERIES FORM OF N(L) 

Substituting the power series (35) for (d/dL) P(V ^ L) of Part 1 and 
the power series (81) for V+(L) into equation (77) yields 

N(L) = tH(-1, 0)a L 2 "- 1 + xff(-l, 0)o 1 L 2 " + • • • (85) 

In the deep-fade region, the leading term dominates the power series 
(85). Therefore, 

iV(L)^T#(-l,0)a L a ' 1 - 1 , for small L; (86) 

■■• N(L) cc L 2 "" 1 , for small L. (87) 

Equations (85) and (86) show that if a -^ and if /i < 1/2, then 
N(L) — > oo as L — > 0. This seems to be nonphysical. Therefore, we 
require that n ^ 1/2. This is the reason we impose this condition on the 
order of singularity of f(a, /3) in Section IV of Part 1.* 

5.1 Prevailing Power Law of Number of Deep Fades 

For the nondiversity fading signals of most radio links, the probability 
density function /(a, /3) of the resultant interfering signal is a smooth 
function which is neither singular nor zero at (a = — 1 , /3 = 0) . Then 
n = 1 and f(a, 0) = H(a, /3). Equations (85), (86), and (87) under this 
condition become 

N(L) = tt/(-1, 0)a L + tt/(-1, 0)a x L 2 

+ [2a rf 4 + W(- 1, 0)a 2 ]L 2 + • • • (85') 

N(L)^TJ(-l,0)a L, for small L. (86') 

N(L) cc L, for small L. (87') 

It is seen that as long as co > /(— 1, 0) > 0, then the expected number 
N(L) of deep fades always obeys the prevailing power law (86'). With 
reference to Fig. 5, the straight lines corresponding to the power law (87') 
have inverse slopes of 20 dB per decade. 



* If we assume that ao = 0, then the only constraint on n is n > due to the unity 
total probability. However, assuming ao = implies that the time derivative of V 
is always zero at V = 0. Such an assumption is unreasonable for multipath inter- 
ference fading, but may be useful if the effects of random circuit interruptions, such 
as equipment failure, are included in the signal fading problem. 
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Apparently, the simple condition is easily satisfied by the nondiversity 
fading signals of most microwave radio links because the experimental 1 _t ' 
data of N(L) are mostly characterized by the inverse slope of 20 dB 
per decade in the deep-fade region. 

As L increases from zero towards unity, equation (85') indicates 
that there may be a transition point beyond which high-order terms 
become significant and the slope begins to deviate. 

For most microwave radio links, the transition points of N(L) seem 
to be well above —10 dB. However, our latest experimental data show 
that the first transition point of N(L) of a relatively short path (path 
length 15.87 miles, / = 4 GHz) is below —20 dB. The inverse slope 
of N(L) in the region from — 10 dB to -20 dB is approximately 20/3 dB 
per decade of number of fades. This indicates that the third term 
[2a d 4 4- 7ra 2 /(-l, 0)]L 3 dominates in the region 0.3 > L > 0.1 for 
this path. 

VI. AVERAGE DURATION OF DEEP FADES 

Substituting the power series (31) of P(V ^ L) of Part 1 and the 
power series (85) of N(L) into equation (79) yields 

iMi=i±m L *» + diL *>+* + d 6 L 2 " +i + . . . 

1{L) = wH(-l, OKI 2 "" 1 + rH(-l, 0)a,L 2 " + ^ 

In the deep-fade region, equation (88) becomes 

l(L)^ — L, for small L; (89) 

fj.a 

.-. t(L) cc L, for small L. (90) 

On Fig. 6, the straight lines corresponding to the power law (90) have 
an inverse slope of 20 dB per decade of fade duration. The experimental 
data 1-6 agree with this conclusion on the slope of l(L) when plotted 
on Fig. 6. 

VII. INVARIANCE OF POWER LAW OF AVERAGE FADE DURATION* 

Equations (36), (87), and (90) show that in general P(V ^ L),N(L), 
and t(L) obey the following set of power laws of deep fades: 

* In equation (26) of Ref. 32, Rice has already predicted that the power law 
t(L) cc L for small L may be applicable to cases more general than the complex 
Gaussian model even though most of his work in Ref. 32 is devoted to the statistics 
of a sine wave plus a narrowband Gaussian noise. 
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P(V ^ L) cc L 2 

N(L) « L 2 

l(L) cc L 



for 



^h 



and small L. 



(91) 
(92) 
(93) 



It is seen that the power laws of P(V ^ L) and N(L) depend on the 
value of n which depends on whether f(a, 0) is smooth, singular, or zero 
at ( a = — 1,0 = 0). On the other hand, the power law (93) for the 
average fade duration is invariant with respect to /z. Since the behavior 
of j{a, 0), and hence the value of n, depends on fading environment, 
we conclude that the power law (93) for the average fade duration is 
insensitive to the fading environment in contrast to the power laws 
of P(V ^ L) and N(L). 

In equation (89), notice that na , and hence t(L), does depend on the 
fading environment. However, it is the power law, t(L) « L, which is 
insensitive to the fading environment. 

Vigants, 4 ' 5,8 Crawford, Hogg, and Kummer 12 have investigated the 
effects of diversity on P(V ^ L), N(L) and t(L). The theoretical results 
and the experimental results of these authors show that in the deep-fade 
region, the diversity drastically changes the power laws of P(V ^ L) 
and N(L) but does not affect the power law t(L) « L. For example, the 
results of Vigants are shown in Table I. From this table it is seen that 
the power laws of P(V ^ L) and N(L) depend on the diversity combina- 
tion of fading signals, whereas the power law, t(L) « L, of average fade 
duration is invariant. 

VIII. INCOMPATIBILITY BETWEEN ONE-ECHO MODEL AND OVERLAND RADIO 
LINKS 

Equations (60), (61), (91), (92), and (93) show that p(L), N(L), and 
1(L) of the one-echo model (with equal magnitudes, A = 1) in the 



Table I — Effects of Diversity on Power Laws 
of Deep Fades 





Nondiversity 


Diversity 


P(V < L) 
N(L) 


cL 
(1/c) L 


d/g) L* 

(2c/ff) L* 

(l/2c) L 



Remark: In this table, the parameter, c, as defined by Vigants, is equal to a 
of this paper; and the parameter, q, as defined by Vigants, is equal to 2/wff( — 1, 0) 
of this paper. 
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deep-fade region are 



u] 



p(L) oc L' 

N(L) cc L° 

1{L) « L J 



for small L. 



(94) 
(95) 
(96) 



On the other hand, the long-term experimental data of a nondiversity 
signal of most overland microwave radio links indicate that 



p(L) « L 

N(L) <x L 

l(L) a L\ 



r for small L. 



(97) 
(98) 
(99) 



The experimental results (97) and (98) disagree with (94) and (95) of 
the one-echo model. 

In view of this disagreement, we may want to check the effect of the 
assumption (81) of V+(L) on the theoretical results of the one-echo 
model. Although we know that the constant term a of (78) generally 
does not vanish, yet we may deliberately set a = and see what kind 
of theoretical results we get. 

If we do so, the theoretical results of the one-echo model become 



P (L) cc L° 




(100) 


N(L) o: L 


r for small L. 


(101) 


l(L) oc L°. 




(102) 



Under this modified assumption, p(L) and t(L) of the one-echo model 
disagree with the experimental results (97) and (99). Similarly, forcing 
the coefficients of other higher order terms of (81) to zero also yields 
theoretical results which disagree with the experimental results. There- 
fore, we conclude that the one-echo model is not suitable for the study 
of the fading signals of most overland microwave radio links. 

However, we emphasize that the experimental data mentioned in this 
section are restricted to the long-term data of overland microwave radio 
links. Therefore, the incompatibility of the one-echo model with these 
data does not necessarily exclude the use of this model for the study of 
other fading problems. 
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IX. APPROXIMATE AVERAGE POSITIVE DERIVATIVE ~V + OF LINE-OF-SIGHT 
MICROWAVE RADIO LINKS 

For line-of-sight microwave radio links, let V tef be the signal level 
when there is no interference, and V, ad (t) be the random fading signal 
when the interference appears. In our analysis 

V(t) = %^ (103) 

' ref 

is a normalized fading signal. 

By comparing the experimental data of t{L) and the theoretical 
equation (89) for t(L), we can estimate the value of a„ of the radio link. 
The value of n in equation (89) can be determined from the experimental 
data on the power laws (36) and (86) for P(V g L) and N(L) of the 
same radio link. 

Our experimental data of several line-of-sight microwave radio links 
in Ohio and Georgia indicate that the value of a ranges from 2 X 10~ 3 
to 4 X 10~ 3 . In the deep-fade region where L is small, equation (81) 
shows that V+(L) = a . 



I™ 



S a = 2 X HT 3 ~ 4 X 10" 3 . (104) 

Substituting (103) into (104) yields 

a* (2 X 10 -3 ~ 4 X 10~ 3 )- V iel . (105) 



it ".-« 



Thus, the average positive derivative of the unnormalized fading signal, 
Vfad(0> of these microwave radio links ranges from 2 X 10~ 3 to 4 X 10~ 3 
times V Tef per second. 

These approximate values of average positive derivative are valid 
only in the deep-fade region because they are deduced from the experi- 
mental data of deep fades. The path length of these radio links ranges 
from 15 miles to 36 miles. The operating frequencies are in 4-GHz 
and 6-GHz bands. 

X. GENERALIZED ASSUMPTION ON V + (L) 

In Part 1, we indicated that the assumption of the Taylor series 
expansion of H(a, /3) is not strictly necessary for the validity of the 
power law of P(V ^ L) deduced from the experimental data. In this 
section, we point out that the assumption of the Taylor series expansion 
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of V+ (L) in Section IV of this part is also not strictly necessary for the 
validity of the power laws of N(L) and t(L). From a theoretical view- 
point, the assumption of f + (L) can be generalized to the following form: 

f + (L) = a + 0(L"), t? > (106) 

where 0(1/') is a symbol to denote the component which goes to zero 
at a rate equal to or faster than that of V as L — * 0. 

In assumption (106), we do not require the existence of the limits 
in equations (83), (84), etc. Therefore, the assumption (106) is less 
restrictive than the assumption (81). It can be shown that the power 
laws of N(L) and t(L) of deep fades based on (106) are the same as 
those based on (81). However, at the present time, we do not have any 
practical evidence to necessitate the use of (106). Therefore, we merely 
point out the possibility but do not explicitly carry out this generalized 
analysis. 
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Part 3. Special Topics on Statistics of a Fading Signal 
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I. INTRODUCTION AND SUMMARY 

In Parts 1 and 2, the analysis was oriented towards an explanation 
of the experimentally observed common behavior of a fading signal TV*. 
The basic assumptions of the theoretical model are kept to a minimum 
in order to include the widest possible variation in practical fading 
environments. During the development of this general analysis, we have 
gained a new insight into several topics related to fading signals as 
investigated by previous authors. 
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Part 3 of this paper is a collection of theoretical treatments of several 
special topics relating our generalized analysis to the work of previous 
authors. These topics include the sum of n unit vectors with random 
phases; a sine wave plus Gaussian noise, ^-distributions, chi-distribu- 
tion, Rayleigh distribution, and log normal distribution. 

In summary: 

(i) In Section II of this part, the results of Part 1 are applied to 
find the amplitude distribution of the sum of n unit vectors with 
uniformly distributed random phases. For n ^ 3, the analysis 
shows that the amplitude distribution always follows the square 
law P(V ^ L) « L 2 for small L. On the other hand, when 
n = 2, the amplitude distribution follows the power law 
P(V gl)«l for small L. 

(it) In Section III, we investigate the model of a sine wave plus a 
narrowband Gaussian noise for the fading signal. By using the 
closed-form solutions of Rice, it is shown that if the power 
spectrum of the Gaussian noise is symmetric with respect to 
the frequency of the sine wave, then the amplitude V and its 
time derivative V are independent; and the conditional average 
positive derivative V+(L) is a constant. On the other hand, if 
the power spectrum is not symmetric, then V and V are de- 
pendent; and the conditional average positive derivative V + (L) 
is a function of signal level V = L. 

As an example, the fading signal spectral density of a mobile 
radio is generally not symmetric with respect to the received 
carrier frequency. Therefore, in the analysis of N(L) and l(L), 
it is not safe to assume that V and V are always independent. 

(in) In Section IV, we investigate the theoretical condition (147) 
on the joint probability density function /(a, /3) of the interfering 
vector such that the amplitude distribution of the fading signal 
belongs to the family of ??i-distributions which includes normal 
distribution, Rayleigh distribution, Maxwell distribution, and 
all of chi-distributions as special cases. 

It is also shown that the set of m-distributions behave like a 
log normal distribution within a small range (148) of signal 
level near its rms value. This result shows that in the interpre- 
tation of the experimental data, one must be cautious in attempt- 
ing to estimate the tails of the distribution by an extension from 
the middle section of the distribution. 

(iv) We find that in general, the integral transformation (20) from 
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j(a, /3) into P(V ^ L) is not unique. Physically this means the 
signals of fading environments with different f(a, /3) can have 
the same amplitude distribution P(V ^ L). As an example, this 
nonuniqueness shows that specifying a Rayleigh distribution 
for P(V ^ L) does not necessarily imply that there are a large 
number of interfering signals; nor does it necessarily imply that 
the real and imaginary parts of the fading signal are normally 
distributed with zero mean. 

II. SUM OF 11 UNIT VECTORS WITH UNIFORM RANDOM PHASES 

The amplitude distribution of the sum of two unit vectors with 
uniformly distributed random relative phase has been shown in Sec- 
tion IX of Part 1 to be, for g L ^ 2, 

P(V gL)=- cos" 1 (l - t) - < 107 > 



7T 

and 

p(L) = ~P(V ^L) = \ 



•[-(-«■. 



(108) 



In the deep-fade region where L is small, this amplitude distribution 
obeys the power law 

P(V £ L) « L. (109) 

The sum of n unit vectors with uniformly distributed random phases 
has been investigated previously by many authors. 33 ,34 ~ 37 The mathe- 
matics involved in obtaining the amplitude distribution for any arbitrary 
n ^ 3 is fairly complicated. Computer numerical integration is needed 
to show the distribution explicitly. In this section, we shall avoid the 
complicated mathematics and shall apply the results of Part 1 to show 
that the amplitude distribution for any arbitrary n ^ 3 in the deep-fade 
region always follows the square law: 

P(V ^ L) <* L 2 , for small L. (110) 

The sum of n unit vectors with random phases can be written as 

Ve'* = jy* (111) 

I = 1 

= [i + |y<<-- 8 *> J '•« (H2) 
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= [l+Re'W 1 , (113) 

where 

Re'' = g «*«'-«. (114) 

i=2 

Rosenbaum 38 has indicated that if all the phases {0,}-:? of the unit 
vectors are independently and uniformly distributed in (0, 2ir), then 
Ve 1 * and Re" have circular symmetric probability density functions; i.e., 
the amplitude and the phase are independent and the phase is uniformly 
distributed in (0, 2tt). 

It then follows that the random signal represented by equation (113) 
for any arbitrary n ^ 3 is a special case of Appendix B. The case for 
n = 2 is an exception because the joint probability density function 
q(R, d) contains a delta function whereas the q(R, 0) in Appendix B is 
assumed to be a smooth function. 

The sum of three unit vectors can be considered as a unit vector 
suffering interference by a random vector R(t)e' Hn which is the sum of 
the other two unit vectors. The amplitude distribution g(R) of R(t) is 
given by equation (108) except for the replacement of the notation L 
by R. Equation (108) implies 

„m = 4- (us) 

Therefore, g(R) for this case is a smooth function which is neither 
singular nor zero at R = 1. Then equation (165) shows that 

P(V ^ L) ~ foCDIr 1 - ^ V « L\ (116) 

The sum of n unit vectors, n ^ 3, can be considered as a unit vector 
suffering interference for a random vector R(t)e iBU) which is the sum of 
the other (n - 1) unit vectors. It is obvious that g(l) ^ simply 
because each of the (n - 1) unit vectors has unity amplitude. Then, 
the results (see Appendix B) imply 

P(V g L) « L 2 , for small L and n ^ 3. (117) 

On a log-versus-dB graph paper, as shown in Fig. 1, the power law (109) 
implies a straight line with an inverse slope of 20 dB per decade of 
probability whereas the square law (117) implies a straight line with 
the same inverse slope of 10 dB per decade of probability as that of 
Rayleigh distribution. 

Therefore, we conclude that for n = 2, the distribution of deep fades 
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is characterized by the inverse slope of 20 dB per decade of probability 
whereas for any n ^ 3, the distribution of deep fades is always charac- 
terized by the inverse slope of 10 dB per decade of probability. This 
conclusion agrees with the numerical results of Norton, et al., in Fig. 2 
of Reference 33. 

III. A SINE WAVE PLUS A GAUSSIAN NOISE 

The statistical behavior of a sine wave plus narrowband Gaussian 
random noise has been investigated in great detail by Rice. 26-27,32 In this 
section we shall apply our general analysis to this case to show the 
consistency of our results with the work of Rice. Furthermore, we shall 
also use the closed-form solution of N(L) and t{L) obtained by Rice 
to show that the conditional average positive derivative V+{L) can be 
either a constant or a function of signal level L, depending on whether 
the power spectrum of the noise is symmetric or asymmetric with 
respect to the frequency of the sine wave. 

In this model, the interfering vector, Re' 9 = a + j/3, represents the 
envelope of a narrowband Gaussian noise; the constant vector represents 
the sine wave with a constant amplitude and frequency /, . The joint 
probability density function f(a, /3) is a two-dimensional normal density 
function; i.e., 

/(«, 0) = 7TT exp [-(a 2 + /3 2 )/26 ] , (118) 

where a and /3 are assumed to be independent normal random variables 
with the same variance b„ and zero mean. 
The well known Rice distribution for the amplitude of this model is 

P ( i ) = ^P(K S /.)=|/.(f)e X p(=^), (.19) 

where 7„(~) is the modified Bessel function of zero th order, and Q is 
the magnitude of the sine wave. In our analysis, Q = 1 because all the 
signals are normalized to the magnitude of the constant vector. 

Rice 27 has also shown that the joint probability density function 
p(Y) V) f° r this model is 

p(t ' v) = ^vm L exp tes mr ~ 2VQ cos * + Q2) 



+ (KV + b&smtfnd*, (120) 
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where <f> is the phase of the resultant fading signal Ve'*; /„ is the frequency 
of the sine wave, w{j) is the power spectrum of the Gaussian noise, and 



k = (271-r P urn - 

Jo 


/,)" df, 


n = 0,l, 2, 


(121) 


B = b b 2 — &i • 






(122) 



3.1 Amplitude Distribution of Deep Fades 

The normal density function (118) is obviously a smooth function 
which is neither singular nor zero at the infinite fade point (a = — 1, 
/3 = 0). Then the results of Section VI of Part 1 predict that the ampli- 
tude distribution of the fading signal in the deep-fade region is 

P(V g L) S tt/(-1, 0)L 2 (123) 

=wM^} L '- (i24) 

On the other hand, the limiting form of the Rice distribution for small L 

is 

■« a W*F)' (i25) 

•'• p(F - L) -2i exp (i6)' L2 ' for l «q = i - ( i26 > 

It is seen that our result (124) agrees with the Rice distribution in the 
deep-fade region. 

The square law (124) implies that on a log-versus-dB graph paper, 
the Rice distribution in the deep-fade region is always characterized 
by the prevailing inverse slope of 10 dB per decade of probability. The 
numerical results of Norton, et al., in Fig. 5 of Reference 33 agree 
with this prediction. 

3.2 Symmetric Power Spectrum, and Constant V+(L) 

If the power spectrum iv(f) of the Gaussian noise is symmetric about 
/, , then bi = and the integration of (120) under this condition yields 

- Pxtf I V)p 2 (V) (128) 

= p 3 (V)p 2 (V), (129) 
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where 

is the Rice distribution for V; 

p ^ ) = vk exp K] (131) 



is a normal density function for V; and 

o 2 = f = V 2 rm . when 6, = 0. (132) 

Oo 

Equations (127) and (129) show that if the power spectrum is symmetric 
about /„ , then the envelope V and its time derivative V are independent, 
and V is normally distributed. 
Substituting (131) and (132) into the definition (75) for V + (L) yields 

f + (L) = =^5? - a = constant. (133) 

V2ir 

Therefore, the conditional average positive derivative V+(L) for this 
model is a constant if the power spectrum is symmetric about /, . 

Substituting (130) and (133) into the general expression (76) for 
N(L) yields 

= la p a (L). (135) 

Then 

.... P ( V S L) 2 f^)^ 

m ~ -mET ~ ^Si) — (136) 

In the deep-fade region where L <£ Q = 1, 

"^vSi^Kl 1 (137) 

^7ra /(-l,0)L, (138) 

l(L) S - L. (139) 

On 
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It is seen that equations (138) and (139) agree with equations (86') 
and (89) of Part 2. 

3.3 Asymmetric Power Spectrum and Nonconstant V+(L) 

If the power spectrum w(j) is not symmetric about /„ , then V and V 
are dependent and &i ^ 0. The joint probability density function 
P(V> V) for this case cannot be written as the product of the individual 
probability density functions of V and V. For this case, Rice 27 has 
obtained N(L) by substituting (120) into the general expression (72) 
and carrying out the integration. This gives 

IB 



N(L) = - ^ ^ l_(~b y 



^ r (QL\^T-n\\ QL 






where 7„(~) is the modified Bessel function of order n, and 



2 _ W 



Comparing equation (140) and the general expression (76) for N(L) 
shows that 






It is seen that when V and F are dependent, then the conditional 
average positive derivative V+(L) is a function of signal level V = L. 
The expected number of fades N(L) and the average fade duration 
t(L) for this case in the deep-fade region are 

iV(L)^7ra /(-l,0)L 
l(L) S £ L, 
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where 

a = lim 7 + (L) = {^J 2, ^j* L. 1 + 2(n + 1)J" 

The work of Clarke, 29 Ossanna, 30 and Gans 31 on mobile radio indicates 
that the power spectrum of the fading signal is generally not symmetric 
with respect to the received carrier frequency unless the straight line 
joining the base station and the mobile antenna is perpendicular to the 
velocity of the mobile and the antenna pattern is symmetric with 
respect to this straight line. Therefore, in the theoretical work of N(L) 
and t(L), it is not safe to assume that V and V are always independent. 

IV. W-DISTRIBUTIONS, CHI-DISTRIBUTIONS, AND RAYLEIGH DISTRIBUTION 

In the study of the experimental data of amplitude distributions of 
short-term high-frequency long-distance propagations, Nakagami 39 
found that the set of experimental data can well be described by a 
family of ?n-distributions:* 

«»-$& !l! ~-[ = ¥\> (143) 

where fl is the mean square value of the fading signal. The operating 
frequency ranged from 10 MHz to 20 MHz and the path length from 
1500 kilometers to 9000 kilometers. Nakagami indicated that these 
results were obtained from short records of data from three to seven 
minutes in length in order to avoid the effects of slow fading on the 
distribution of rapid fading. 

The various properties of the m -distributions have been investigated 
in detail by Nakagami. 39 It is easily shown that the set of chi-distribu- 
tions 28 is a subset of w-distributions by setting 2m = any positive 
integer in (143). This means the normal distribution, Rayleigh dis- 
tribution, and Maxwell distribution are also special cases of m-distri- 
butions when m = 1/2, 1, and 3/2 respectively. On Rayleigh paper, 
all the m -distributions appear to be straight lines passing through the 
common point of 50 percent at dB, with different slopes which depend 
on the value of m. The graphical representation of ?n-distributions can 
be seen in Reference 39. 

However, in Reference 39, one does not know the theoretical condi- 
tion under which the amplitude distribution of a fading signal will 



* To avoid possible confusion, we emphasize that the Nakagami distribution 
mentioned in Section I of Part 1 is not the m-distribution discussed in this section. 
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belong to this family of ?w-distributions. In this section we shall find 
the condition on the joint probability density function f(a, /3) of the 
interfering vector such that the amplitude distribution P(V ^ L) will 
belong to m-distributions. 

Expanding the exponential function in equation (143) into a power 
series gives 

Comparing equation (144) and the general power series (35) for p(L) 
shows that 

/x = m, (145) 

(28 + 2 M ) d 2S+2 = ,* ( ~ l) " (146) 



Substituting (34) into (146) gives 

JT_ y» H 2S - 2 ,, 2 ,(-l, 0) (-1) a = i o • • • (147) 

Thus, equations (145) and (147) are the general conditions on the 
interfering vector such that P(V ^ L) is an m-distribution. 

In Section V of Part 2 we showed that /x ^ 1/2 which implies m ^ 1/2. 
Nakagami 39 has also found this condition on the parameter m by a 
different approach. Since m can be any value ^ 1/2, equation (143) 
represents an infinitely large family of distributions. 

4.1 Log Normal Behavior of m-Distribution Near the RMS Value 

The experimental data of optical propagation 40-42 and line-of-sight 
radio links show that the distributions of the signal scintillation near 
its average value are approximately log normal. Usually the accuracy 
of experimental data is best in the middle section of the distribution 
and deteriorates towards the tails. It is quite tempting to estimate the 
tails of the distribution by an extension from the middle section. The 
deviations of the experimental data at the tails are often attributed to 
the experimental error. 

However, de Wolf 40 and Deltz and Wright 42 have pointed out that 
the use of the middle section of a log normal paper may not be a reliable 
test of the log normal distribution. The differentiation between the log 
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normal distribution and certain other distributions may be significant 
only at the tails rather than the middle section of these distributions. 
Nakagami 39 has pointed out that all the m-distributions behave like 
a log normal distribution for_the fading signal V(t) in the neighborhood 
of its rms value <r R = Vfi. The explicit bounds on the signal level 
within which this approximation holds are: 



20 



logio 



Vfi 



« 4.3 dB. (148) 



It is shown in Appendix F that the ^^-distribution within the signal 
range (148) is approximately equal to 



P(L) = 



2m' 



_L-T(m) 



exp [-2m(ln L - In Vty 2 ], (149) 



which is a log normal distribution for the signal level L. 

Therefore, the m-distributions, including the normal distribution, 
Rayleigh distribution, Maxwell distribution, and chi-distributions, all 
behave like a log normal distribution within the signal range (148). This 
result points out that in the interpretation of the experimental data, 
one must examine the behavior of the data not only inside but also 
outside of the range (148) in order to assert their distribution. 

V. NONUNIQUE RELATION BETWEEN AMPLITUDE DISTRIBUTION AND f(a, /3) 

From the general integral relation (20) between P(V ^ L) and 
f(a, /3), it is seen that any component of /(a, j3) that is antisymmetric 
with respect to (1 + a) and/or (3 will cancel out in the integration (20), 
and contributes nothing to P(V ^ L). This means that there are many 
different f(a, /3)'s, with the same symmetric* part and different anti- 
symmetric* part, which correspond to the same amplitude distribution 

P(V ^ L). 

Mathematically, this means the integral transformation (20) from 
f(a, /3) into P(V ^ L) is not unique. Physically, this means the fading 
signals in fading environments with different /(a, 0) can have the same 
amplitude distribution. 

Furthermore, even if we restrict /(a, 0) to functions symmetric with 
respect to (I + a) and 0, the relation between j(a, /3) and P(V ^ L) is 
still not unique. We shall demonstrate this nonunique relation specifi- 
cally by using the results for the m-distributions previously discussed. 
We notice that for each S, equation (147) is an algebraic equation 

* With respect to (1 + a) and 0. 
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for (S + 1) unknowns \H 2S -2p,2»(— 1, 0)}"„=o • It is then obvious that 
there are infinitely many different sets of {H 2S -2v,2v(— 1, 0)}"„1o which 
will satisfy equation (147) because there is only one equation for 
(S + 1) unknowns. This nonuniqueness gives a great freedom for the 
wide variations of the individual term, H 2 g-2v,2w(—l, 0), which is the 
even-order partial derivative of H(a, /3) at (a = — 1, j8 = 0). From the 
Taylor series (27), it is seen that H 2S -2y,2»{— 1, 0) is the coefficient of 
the even-order term (1 4* ot) 2S2v {S 2 " which is symmetric with respect 
to (1 4- oc) and /3. Therefore, the relation between /(a, /3) and P(V ^ L) 
is not unique even if f(a, /3) is symmetric with respect to (1 + a) and /3. 
A more detailed discussion of this nonunique relation in polar co- 
ordinates is given in Appendix G. 

VI. PHYSICAL MODEL AND RAYLEIGH DISTRIBUTION 

In this section, we shall show that specifying a Rayleigh distribution 
for the amplitude of a complex fading signal TV* does not necessarily 
imply that there is a large number of interfering signals; nor does it 
necessarily imply that the real part and the imaginary part of the 
fading signal are normal with zero mean. 

Let x and y be the real part and the imaginary part respectively of 
the complex fading signal Ve'*, and let F(x, y) be the joint probability 
density function of x and y. Since V 2 = x 2 + y 2 , then the probability 
of V ^ L is the probability of x and y falling within the circular region 

x 2 + y 2 ^ L 2 . (150) 

Therefore, P(V 5= L) is the integration of F(x, y) over the circular 
region (150); i.e., 



P(V £L)=[" f __^_ F(x, y) dx dy. (151) 

J„=-i, J«=-Vr.*-v* 



6.1 Number of Interfering Signals 

In the most common derivation of the Rayleigh distribution, the 
fading signal is assumed to consist of a large number of random, inde- 
pendent interfering signals, 

Ve" = E Eie 1 " - x + jy. (152) 

Furthermore, it is assumed that none of the components {-E,})!" pre- 
dominates in the summation (152). Then by the central limit theorem, 
one argues that as the number, n, of interfering signals approaches 
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infinity, the real part x and the imaginary part y of the fading signal 
Ve'* become independent normal random variables with the same 
variance and zero mean. This implies that V and <j> are independent and 
<t> is uniformly distributed in (0, 2ir). Under this condition, the distribu- 
tion of the random amplitude V is Rayleigh. 

However, by an observation similar to those in Section V of Part 3 
and Appendix G, we realize that the transformation (151) from F{x, y) 
into P(V ^ L) is not unique. Given an amplitude distribution P(V ^ L), 
there correspond infinitely many different F(x, y)'s. In other words, 
the independent normal distribution for x and y with the same variance 
and zero mean is only a sufficient condition but is not a necessary 
condition for the amplitude V to have a Rayleigh distribution. Since 
x and y do not have to be normal, then the number n of the interfering 
signals does not have to be large. 

In Appendix G, we have shown that the relation between P(V ^ L) 
and F(x, y) becomes unique if the following two additional conditions 
are imposed; 

(i) V and <f> are independent, and 
(ii) <f> is uniformly distributed in (0, 2tt). 

For long radio links such as beyond-the-horizon radio links, the condi- 
tions i and ii seem to be applicable. However, for line-of-sight radio 
links, our experience indicate that the phase <t> has much higher tendency 
of wide variation during the deep fade where V is small. This means 
for short radio links, V and <£ may not be independent and <t> may not 
be uniform. Therefore, in our general analysis we do not impose the 
conditions i and ii. 

6.2 Mean Values of x and y 

In the integral relation (151) the antisymmetric part of F(x, y) 
contributes nothing to the amplitude distribution P(V ^ L), but does 
affect the mean value of x and y. Then by adding a suitable* antisym- 
metric function to F(x, y), the mean values of x and y can be changed 
arbitrarily without affecting the amplitude distribution P(V ^ L) of 
the fading signal Ve'*. 

In other words, given an amplitude distribution P(V fS L) of a 
complex fading signal Ve'*, the mean values and the higher moments 



* The probability density function F(x, y) must be ^0 for any x and y; therefore, 
the symmetric part of F(x, y) must be ^ | antisymmetric part of F(x, y)\ when the 
antisymmetric part is negative. 
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of x and y are not unique. (However, the moments of the amplitude 
V are unique.) 

Therefore, specifying a Rayleigh distribution for the amplitude V 
does not necessarily imply that the mean values of x and y are zero. 
Physically if there is a direct path between the transmitter and the 
receiver of a radio link, then the mean values of .t and y may not be zero. 
However, the results of this section show that the mere nonzero means 
of x and y do not necessarily exclude the Rayleigh distribution for the 
amplitude of the fading signal. 
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APPENDIX A 



List of Symbols and Their Definitions 

A The constant amplitude of the echo in the one-echo 

model. 
a The zero th order term of the Taylor Series expansion 

of V + (L) denned in equation (81). 
a n The coefficient of the rath order term of the Taylor 

series expansion 7+(L) defined in equation (81). 
B = b b 2 — b\ as defined in equation (122). 

b n ,n = 0,1,2 Defined by equation (121). 
C 4 The coefficient of the fourth order term L 4 of P(V ^ L) 

in equations (161) and (163). 
C" Defined by equation (30). 

d» s+ 2 The coefficient of the power series representation of 

P(V g L) defined in equations (31) and (34). 
f(a, /3) Joint probability density function of a and /3. 

F{x, y) Joint probability density function of the real part x 

and the imaginary part y of the fading signal. 
f q The frequency of the sine wave. 

g(R) The probability density function of the amplitude R of 

the interfering vector. 
g n (R) The nth order derivative of g(R). 

H(a, /3) The smooth part of /(a, 0) as defined in equation (23). 

H n , m (a, (3) The partial derivative of H(a, 0) as defined in equa- 

tion (29). 
H(-l, 0) The value of H(a, p) at (a = -1, = 0). 

H(-l, 0) The average value of H(a, 0) at (a = -1, P = 0) if 

H(a, /3) is discontinuous at this point. 
I„(~) Modified Bessel function of order n. 

L An arbitrarily specified signal level in the study of the 

statistics P(V ^ L), N(L) and t(L). 
m An integer. 

n An integer. 

N(L) Expected number of fades per unit time below the 

specified signal level L. 
0(L") A symbol to denote a function which goes to zero at a 

rate equal to or faster than V as L — > where 77 > 0. 
P(V ^ L) Probability that the amplitude V of a fading signal 

fades below a specified signal level L. 
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p(L) Probability density of the amplitude V at 1fhe specified 

signal level L. 

■p 2 (V) Probability density function of V. p 2 (V) = p(V). 

p(V> V) Joint probability density function of V and V. 

j)i(V | L) Conditional probability density of V under the condi- 

tion V = L. 

Q The constant magnitude of the sine wave in the 

specialized model of a sine wave plus a Gaussian noise. 
In this paper, Q = 1 because of the normalization of 
the signal level. 

q(R, 9) Joint probability density function of R and 0. 

R The amplitude of the resultant interfering vector Re' 9 

of the fading signal model. 

Re'" The resultant complex interfering vector. 

r An integer. 

S An integer. 

t A variable representing time. 

t(L) Average duration of fades below L. 

V The amplitude of the envelope of a complex fading 
signal normalized to the nonfaded value V rof . 

V The time derivative of the normalized amplitude V of 
the fading signal. 

Ve'* The envelope of the fading signal. 

V rmK The rms value of the time derivative of V. 

V fad (t) The unnormalized amplitude of the random fading 

signal. 
V tB , The nonfaded signal level when there is no interference. 

V+ (L) Conditional average positive derivative of V as defined 

in equation (75). 
W{6) Probability density function of the random relative 

phase 9 of one-echo model. 
W n {9) The nth order derivative of the probability density 

function W(6). 
x The real part of the fading signal Ve'*. 

y The imaginary part of the fading signal Ve'*. 

a The real part of the interfering vector Re'". 

/3 The imaginary part of the interfering vector Re 1 . 

r(~) Gamma function. 

y 2 Defined by equation (141). 

7] An arbitrary constant > 0. 

9 The phase of the resultant interfering vector Re' . 
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d L Defined by equation (53). 

n (1 — n)/2 is the order of singularity of the joint prob- 

ability density function j(a, 0) at (a = — 1, = 0) as 
defined in equation (23). 

v An integer. 

p (V, 0) Joint probability density function of the amplitude V 

and the phase of the fading signal. 

<t(V, 0) A function of V and satisfying the homogeneous 

integral equation (207). 

The phase of the envelope of the complex fading signal. 

J2 The mean square value of an m-distributed random 

variable. 

APPENDIX B 

Amplitude Distribution In Polar Coordinates 

Since the results of P(V g L) for m = 1 cover a large class of fading 
problems and since the statistical behavior of the interfering vector is 
sometimes described by the joint probability density function q{R, 0) 
of the interfering vector, we shall also obtain the power series repre- 
sentation of P(V ^ L) in terms of q(R, 0) when n = 1. By using the 
relations: 

a = R cos 0, (153) 

= R sin 0, (154) 

and the Jacobian relation "between f(a, 0) and q(R, 0) one can represent 
the coefficients {d M+a j of (42) in terms of q(R, 0). This gives 

P(V ^ L) = rq(l, k)L 2 + d 4 V -f d, L 6 + • • ■ (155) 

where 
d 2 = *q(l } ir) = *■/(- 1,0), (156) 

d* = I [«(1, t) - ff ll0 (l, t) + fc.oU, ») + ffo..(l, t)] etc., (157) 

o 

fc-<fc«-jFSF «*'«>■ (158) 

In the deep-fade region, 

P(V ^ L) > 7r 5 (l, tt)L 2 , (159) 
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and 

p(L) > 2irq(l, r)L. (160) 



i-o 



B.i Circular Symmetric Probability Density Function 

In this subsection, we consider a special case where the interfering 
random vector, Re' e , has a circular symmetric probability densitj' 
function q(R, 0); i.e., R and 8 are independent and 6 is uniformly dis- 
tributed in (0, 2t). For this case, let g(R) be the probability densitj' 
function of the magnitude R of the interfering vector, then 

P(V ^ L) = i ff (i)L 2 + CJJ + CJ? + • • • , (161) 

q(R, 6) = g(R) ^ , (162) 

Ci = w [9{1) ~ 9l(1) + 92{1)] etc -' (163) 

9n(D = ^ g(R) U-, • (164) 



In the deep-fade region 

P(V SL) >lgO)L\ (165) 

L->Q 

These results are used in Part 3 where we discuss the relation between 
our generalized analysis and the existing theoretical work. 

appendix c 

Nonanalytic H(a, 0) 

In Section V of Part 1 , the analysis is carried out based on the assump- 
tion that H(a, 0) can be expanded into a two-dimensional Taylor series. 
In this appendix, we shall investigate two cases where H(a, /3) cannot 
be expanded into the two-dimensional Taylor series. The objective is 
to show that from a theoretical viewpoint, the assumption of the Taylor 
series expansion of H(a, /3) is not strictly necessary for the validity of 
the power laws of deep fades discussed in this paper. 

c.i Continuous H{a, 0) With Unbounded Derivatives 

If H(a, 0) is continuous at (a = — 1, = 0) but its first-order partial 
derivatives and/or its higher order partial derivatives are unbounded 
at (a = — 1, /3 = 0), then H(a, &) in the neighborhood of a = —1 and 
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/3 = can be written as 

H(a, 0) = H(-l, 0) + 0{[(1 + a) 2 + PY"}, (166) 

where 1 > 77 > 0, and 

0|[(l+a) 2 + /3T /2 } 
is a symbol to denote the component which goes to zero at a rate equal 
to or faster than that of 

[(1 + a) 2 + &Y 2 as [(1 + a) 2 + /3 2 ] -> 0. 

It is obvious that H(a, 0) given by (166) cannot be expanded into 
Taylor series because the derivatives of H(a, 0) are unbounded (i.e., 
singular) at (a = — 1 , = 0) . 

Substituting equations (166) and (23) into the general formulation 
(20) for P(V ^ L), and carrying out the integration yields 

P(V ^ L) = 7r/f( ~ 1, ° ) L 2 " + 0(L 2l,+ ") t (167) 

p(L) = 2nH(-l, OL 2 "" 1 + 0(L 2 *- 1+ '), (168) 

where 0(L 2 " + ") is a symbol to denote the high-order terms which go to 
zero at n rate equal to or faster than that of L 2 " + ' as L — » 0. Since tj > 0, 
then in the deep-fade region 

P( V < L) -^ 7rg( ~ 1 '° ) L 2 " , (169) 

which is the same as equation (36). Then the discussions and conclusions 
in Sections VI, VII, and VIII of Part 1 on the power laws of deep fades 
for n — 1, 1 > /x ^ 1/2 and m > 1 are readily applicable to the present 
case even though the derivatives of H(a, 0) are unbounded at (a = —1, 
= 0). 

c.2 Discontinuous H(a, 0) 

Suppose H(a, 0) and its derivatives are bounded but are discontinuous 
at /3 - so that 

lim //„.„,(«, 0) ^ lim H n . m (a, 0); (170) 

0-(l + 0-0- 

i.e., 

H (a, + ) 5* //„,,„(«, 0") 

n = 0,1,2, ••• , m = 0,1,2, ••• . (171) 
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Then on each side of |S = 0, the one-sided Taylor series expansion of 
H (a, /3) is applicable. One for /3 > and another for < 0. Substituting 
these two Taylor series and equation (23) into the general formulation 
(20) for P(V ^ L) and carrying out the integration, one can show that 

P(V ^ L) > *•#(-! >0) L 2^ (172) 

p(L) > 2tH(-1, 0)L 2 "-', (173) 

Z.-0 

where 

ff(-l, 0) = Hff(-1, + ) + tf(-l, 0-)] (174) 

is the average value of the discontinuous H(a, 0) at (a = — 1,0 = 0). 
It is seen that equation (172) is also the same as equation (36) except 
for the proper interpretation of ff(— 1, 0) when H (a, 0) is discontinuous 
at (a = — 1 , /3 = 0) . Therefore, the discussions and conclusions of 
Sections VI, VII, and VIII of Part 1 are also applicable to the present 
case. 

APPENDIX D 

Integration jor Power Series of Amplitude Distribution 
Substituting (27) and (23) into (20) gives 



P(V g L) = £ -Jj E CT„- r . r (-l, 0)/„_ r . r , (175) 

n = 0'il r=0 



where 



\-L i.o— l+Vz'-fl 



/.-,,, = f [ [(l+ay + ffi-^l+aypdadp. (176) 

Jfi x, J Q __i_Vi»-/3> 

From (176) it is seen that if either (n — r) or r is an odd integer, then 
I n - TT vanishes because the integrand is antisymmetric. Therefore, 
I„- r ,r does not vanish only when both (n — r) and r are even integers. 
Then let n = 2S and r = 2v. Equation (176) becomes 



= f f " — [(l+a) 2 + {S 2 Y-\l+ay s - 2 'p 2 'd*d0 

JB--I, ■'n-l-v'l'-S' 



_ 2TQS - v + flrfr + \) 2S+2M 

(25 + 2n)Y(S + 1) ^ u " ; 
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Combining (175) and (177) gives 
P(V ^ L) = E L 2S+2m 

S-0 

r 2 ^ cifns - v + w(y + h) J7 /i m l 

' |_(2S)! h C2S + 2/*)r(S + 1) rt2S - 2 - 2 ^ x ' u; J 

(178) 
= E<2 2S+2 L 2S+2 ", (179) 



where 



, 2 gg C 2 2 fr(g - y + \)V(y + j) „ , - m 

rf 2S+2 - (2jS) , Z. (2(S + 2/x)r(<s + 1} fli«-i,. a .l i, u;. 



(180) 



APPENDIX E 



Derivation for Expected Numbers oj Fades 

Suppose T = t 2 — U is the time interval in which we want to find the 
expected number of times that the random fluctuating signal V(t) 
crosses the signal level V = L. This interval is divided into a large 
number of smaller intervals of width At so short that each contains 
no more than one level crossing. We first consider the expected number 
of upward level crossings. The downward level crossings can be treated 
similarly. 

In an infinitesimal interval At, the conditions for an upward level 
crossing of V(t) are 

v(o = mp- > o (18D 

at 

V(t) At > [L - V(t)} > (182) 

These two conditions are shown graphically in Fig. 8. On a V versus 
V plane, the region in which V and V satisfy conditions (181) and (182) 
are shown as the shaded area in Fig. 9. The integration of the joint 
probability density p(V, V) over this range will give the probability 
that V(t) will have an upward level crossing in At, 

PUL) - T ° f ■ P(*i V)dVdV (183) 

J K = Jr-L-VAt 



£< At f] " Vp(V, V) 



dV. (184) 

V = L 
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V(t) 

"-— v(D 

\<. LEVEL L 



L-V(t) 
Wt) 



1 



[L-V(t)]>oL^ 




^ At *| 



t t+At 

Fig. 8 — The conditions for an upward level crossing of V(t) in an interval At. 



The expected number N aD (L) of upward level crossings per unit time is 
PUD 



N W (L) = -^ = f Vp(V, V) 



dV. 



(185) 



Similarly, the expected number of downward level crossings per unit 
time is 



N, 



„. wn (L) = f" \V\p(V,V) 

J V--QO 



dV. 



The total expected number of level crossings per unit time is 
NXD = iV up (L) + N dowa (L) 



= f. 1 I V | p(V, V) 



dV. 



(186) 



(187) 




Fig. 9 — The region in which V and V satisfy the conditions for an upward level 
crossing. 
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The expected number of fades per unit time below V = L is 
N (L) = $N C (L) = N UP (L) = N dovn (L) 



= f V ~~ Vp(V, V) 



dV. 



V = L 



(188) 



APPENDIX F 

Log Normal Behavior of M -Distributions Near the RMS Value 
Let 

Z = 20 log 10 (7^=) = 20[log 10 L - log 10 Vn] (189) 

be the signal level in dB with respect to its rms value \/H. Equation 
(189) implies 

(190) 

(191) 

(192) 

(193) 



Z . / L 



= exp 



L 

dL Vtt 



dZ 



M 



exp 



M 



L_ 

M ' 



where 



M = 20 log 10 e = 8.686 dB. 



Then the probability density functions of L and Z are related by the 
Jacobian relation: 28 ' 18 



h(Z) = p(L) 



dL 
dZ 



= ±P(D. 



(194) 



Substituting the m -distribution (143) and equation (191) into (194) 
yields the following probability density of Z. 



h(Z) = 



2m, 



[=¥\ 



M 2m 

MT(m) K^a) GXP 
2m m / \2Z (2Z 

Mf5o exp riM - exp \ M 



Substituting the following power series 

CO 



2Z\ -A 1 2ZY 
exp \—f) = 1^Z\ lT7 



(195) 



(196) 
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into (195) gives 

h(z) = bw~ m J exp r m [ 2 (i) + Sii 



2ZY 

M, 



If the signal level L is very close to its rms value "\/ft> then 



L 



'2Z 
M 



Vtt 



= 2 



= - 1 



«1, 



In 



L 



Va 



«i, 



for 



(197) 

(198) 
(199) 



(200) 



| Z | « M/2 
= 4.3 dB 
Under this condition, we have 

Then h(Z) in equation (197) becomes 

A ^ = [^R e ""] ex p[- 2m (l) 2 ]- (201) 

which is a normal distribution for Z. Substituting (190) and (201) into 
(194) yields 



P(D = 



f 2m m _ m ] 
lLT(m) e J 



exp [-2ra(ln L - In Vo) 2 ], 



(202) 



which is a log normal distribution for L. Therefore, within the range 
| z | « M/2 = 4.3 dB, all the w-distributions for any m ^ 1/2 behave 
like a log normal distribution given by (202). 



APPENDIX G 

Nonunique Relation Between P(V ^ L) and /(a, /3) 

In the definition of the fading signal model discussed in Section II 
of Part 1, the four random variables V, <j>, a, and /3 are related by 

1 + a = V cos <t> J 
/3 = Fsin^j 



(203) 
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Let p(V, (p) be the joint probability density function of the amplitude 
V and the phase <p of the fading signal. By the Jacobian 28 ' 1 of the 
transformation (203), it is easily shown that 

p (7 ) ) = V(l + a) 2 + /3 2 /(a, 0) = Vj(V cos0, 7sin0). (204) 

In this appendix we shall show that the relation between P(V ^ L) 
and p(V, <f>) is not unique. This then implies that the relation between 
P(V ^ L) and /(a, /3) is also nonunique because of the simple algebraic 
relation (204) between p(V, <p) and /(a, 0). 

In the study of the amplitude distribution P(V ^ L) of a complex 
fading signal IV*, it is often assumed that V and 4> are independent 
with <p uniformly distributed in (0, 2?r). However, in the study of 
interference, distortion, FM radio system, radio navigation system, etc., 
many authors" 7 ■■ ,3>44il ° have investigated the distribution of the random 
phase, <p(t). These results show that the distribution of phase is not 
always uniform. Furthermore, when the signal is weak (i.e., V is small), 
the phase is more likely to vary over wider range. This means the random 
variables V and <p are somewhat correlated. Therefore, V and <p generally 
can be either dependent or independent and (p can be either uniformly or 
nonuniformly distributed. 

By definition, 28,18 the probability density p(V) of V is the integration 
of p(V, <p) over the entire range of 0; i.e., 

p(V) = [ * f>(V,<p)d<t>. (205) 

Furthermore, the cumulative amplitude distribution P(V ^ L) is the 
integration of p(V) from V = to V = L. Therefore, 



P(V £L)= [ f P (V,<p)d<pdV. 

JV"0 •'0=0 



(206) 



Given a joint probability density function p(V, <p), then P(V g L) can 
be calculated by (206). 

On the other hand, given an amplitude distribution P(V ^ L), 
equation (206) is an integral equation to solve for p(V, <p). An immediate 
question arising in solving the integral equation (206) is the uniqueness 
of the solution. A procedure to test the uniqueness of the solution is to 
consider the following homogeneous equation 



= f " f T a(V,<p)d<pdV. (207) 

If the homogeneous equation (207) has a nontrivial solution, then the 
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solution of (206) is not unique because, given any particular solution 
Pp (V, 0) of (206), then 

p(V, 0) = p P (V, 0) + ca(V, 0) (208) 

is also a solution of (206) where c is an arbitrary constant.* 
It is obvious that all of the following functions 

*n(V, 0) = UV) sin (n0) + UV) cos (n<t>); n = ±1, ±2, ±3, • • • 

(209) 

are nontrivial solutions of the homogeneous equation (207) where 
£„(F) and f„(F) are arbitrary functions of V. Notice that the nontrivial 
solutions (209) contain both symmetric and antisymmetric functions 
of 0. Furthermore, any arbitrary linear combination of \a n (V, 0)} is 
also a solution of the homogeneous equation (207). By the experience 
of Fourier series synthesis technique, we know that the linear combina- 
tion of the set \<r n (V, 0) { is able to represent a very large class of either 
simple or complicated functions of V and 0. 

Therefore, given an amplitude distribution P(V ^ L), the integral 
equation (206) lias infinitely many different solutions. 

On the other hand, in equation (206), if one imposes the following 
two additional conditions: 

(i) V and are independent, and 
(ii) is uniformly distributed in (0, 2ir), 



^*' = >» = M P( ^ I) 



then 

i 1 a 

(210) 

T UU t = K 

is the only possible solution of (206). 
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